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HIGH LEVELS OF CIRCULARLY POLARIZED EMISSION FROM THE RADIO JET IN NGC 1275 (3C 84) 
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ABSTRACT 

We present multi-frequency, high resolution VLBA circular polarization images of the radio source 3C 84 in 
the center of NGC 1275. Our images reveal a complex distribution of circular polarization in the inner parsec of 
the radio jet, with local levels exceeding 3% polarization, the highest yet detected with VLBI techniques. The 
circular polarization changes sign along the jet, making 3C 84 also the first radio jet to show both signs of circular 
polarization simultaneously. The spectrum and changing sign of the circular polarization indicate that it is unlikely 
to be purely intrinsic to the emitted synchrotron radiation. The Faraday conversion process makes a significant 
and perhaps dominant contribution to the circular polarization, and the observed spectrum suggests the conversion 
process is near saturation. The sign change in the circular polarization along the jet may result from this saturation 
or may be due to a change in magnetic field order after an apparent bend in the jet. From the small spatial scales 
probed here, ~ 0.15 pc, and the comparably high levels of circular polarization inferred for the intra-day variable 
source PKS 1519-273, we suggest a connection between small spatial scales and efficient production of circular 
polarization. 

Subject headings: galaxies : active — galaxies; jets — galaxies: individual: 3C 84, NGC 1275 — polarization 
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(VLBA) circular polarization images of the radio source 3C 84 
in the center of NGC 1275 (z = 0.017). Because of its proxim- 
ity, we have very high spatial resolution, ~ 0.15 pc, in the radio 
jet of 3C 84. We report very strong circularly polarized emis- 
sion in an optically thin region, approximately 0.5 pc south of 
the peak radio emission. In §2 we describe our observations 
and calibration techniques. Our results are presented in §3, and 
we analyze and discuss them in §4. 

2. OBSERVATIONS 

We observed 3C84 in December of 1997 (epoch 1997.94) 
with the VLBA at four frequencies: 5, 8, 15, and 22 GHz. Our 
nearly full track observations were divided between the observ- 
ing bands, giving excellent (u,v)-plane coverage and approxi- 
mately 2 hours on-source per frequency. Here we present and 
discuss the CP results from our 15 and 22 GHz bands. Total 
intensity (Stokes /) images are presented in figure 1 . 

To calibrate the antenna gains to detect circular polarization, 
we used a modification of the Zero-V self-calibration technique 
described by Homan & Wardle (1999). Zero-V self-calibration 
assumes there is no real circular polarization in the data and cal- 
ibrates the baseline correlations {RR = I+V,LL = I-V) versus a 
pure Stokes / model. The derived antenna gains will act to sup- 
press any real Stokes V in the data; however, except in the case 
of a point source, the gain "corrections" cannot remove the CP 
signal completely because the antenna gains are multiplicative 
while V is additive in the baseline correlations. 

For an extended source, like 3C 84, the effect of Zem-V self- 
calibration on real Stokes V will be to reduce its level and in- 
duce V of the opposite sign on other strong structure in the 
source. It is possible to reconstruct the original circular polar- 
ization distribution through careful modeling procedures, and 
that is what we have done. For the 15 and 22 GHz results pre- 
sented here, the core region is well resolved with a slowly 
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1. INTRODUCTION 

Circular polarization from synchrotron emitting radio jets 
can be a powerful diagnostic of the relativistic fluid transported 
by the jets. In principle, circular polarization (CP) observa- 
tions can be used in concert with linear polarization and spectral 
information to constrain the low energy end of the relativistic 
particle distribution, measure magnetic field order and strength, 
and perhaps even deduce the particle content, e^e~ vs. p'^e~, of 
jets (Wardle et al. 1998). 

For inhomogeneous synchrotron radiation, circular polariza- 
tion may be produced as an intrinsic part of the emitted radi- 
ation or through the process of Faraday conversion of linear 
polarization to circular (Jones & O'Dell 1977). On theoreti- 
cal grounds (e.g. Jones 1988) the Faraday conversion process 
is expected to dominate over the intrinsic CP of synchrotron 
radiation in typical jets. 

In recent years circular polarization has been detected in a 
wide range of synchrotron sources: a large number of powerful 
active galaxies (Wardle et al. 1998; Homan & Wardle 1999; 
Rayner, Norris, & Sault 2000; Homan, Attridge, & Wardle 
2001), low-luminosity AGN such as Sagittarius A* and M81* 
(Bower, Falcke, & Backer 1999; Sault & Macquart 1999; Brun- 
thaler et al. 2001; Bower, Falcke, & Mellon 2002), intra-day 
variable sources (Macquart et al. 2000), and galactic micro- 
quasars (Fender et al. 2000, 2002). However, the spectral evi- 
dence from those objects with constraints at multiple frequen- 
cies creates a confusing picture which could be consistent with 
either intrinsic CP or Faraday conversion (Bower et al. 1999; 
Sault & Macquait 1999; Macquart et al. 2000; Fender et al. 
2000; Brunthaler et al. 2001; Fender et al. 2002). Only in the 
case of 3C279 (Wardle et al. 1998) is the Faraday conversion 
process strongly favored by direct observational evidence. 

Here we present multi-frequency Very Long Baseline Array'* 
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distortion will occur from our initial assumption of V = 0. 
3. RESULTS 

Our circular polarization images of 3 C 84 are displayed in 
figures 2 and 3, and we find very strong CP at both 15 and 
22 GHz. The registration between the two frequencies was 
performed on the Stokes / images alone and is illustrated for 
Stokes y in figure 3. Both the positive and the negative features 
appear at the same locations at the two frequencies frequencies 
to within a small fraction of a beam-width. All the circularly 
polarized emission is located south of the map peak, with the 
change in sign from positive to negative occurring after the jet 
appears to bend to the east. 

The bright knot of positive circular polarization (component 
"A" in figure 3) is unresolved at 15 GHz and only barely re- 
solved along the jet at 22 GHz. It is located in a part of the 
jet that has a smoothly declining brightness profile with no dis- 
tinct features in total intensity. This region of the jet is opti- 
cally thin with a spectral index of a = -1 to -1 .5 (5 cx i^^") for 
Stokes / between 15 and 22 GHz. The fractional local circular 
polarization for this spot is Wc = +3.2 ± 0.1% at 15 GHz and 
nic = +2.3 ± 0.2% at 22 GHz. The values are taken at the loca- 
tion of the Stokes V peak, and we can find the spectral index 
of the fractional circular polarization: oc j/-0 9±0 3 ^ote that 
the data were weighted and tapered at the two frequencies to 
give an effectively matched resolution for these measurements. 

The negatively polarized components ("B" and "C" in fig- 
ure 3) are further upstream where the jet is optically thick with 
a ~ +1 .2 at the location of "B" and a ~ +1 .5 at the location of 
"C". These components have fractional circular polarizations 
of nic = -0.7 ± 0.1% at 15 GHz and = -1.3 ± 0.2% at 22 
GHz for component "B" and = -0.6 ±0.1% at 15 GHz and 
Wc = -1 .0 ± 0.2% at 22 GHz for component "C". These values 
include a +10% amphtude correction at 15 GHz and a +30% 
amplitude correction at 22 GHz to compensate for the typi- 
cal reductions introduced by our calibration procedure, see §2. 
The spectral indicies of fractional polarization for these nega- 
tive components are nic for "B" and oc 4±o.7 
for "C". 

It is important to note that our 5 and 8 GHz images (Homan 
and Wardle, in prep.) show no significant circular polarization 
in this core region, although a predominantly negative sign of 
CP is seen farther south in the jet at 5 GHz. The core region 
is poorly resolved at these longer wavelengths and has become 
optically thick, so both blending and optical depth may play a 
role in the lack of detected CP. At the approximate location of 
component "A", the Stokes-1 spectrum is flat between 8 and 15 
GHz and becomes inverted between 5 and 8 GHz, suggesting 
that the self-absorption turnover for this part of the jet occurs 
somewhere between 8 and 15 GHz. Assuming t = 1 at 10 GHz, 
we estimate the optical depths at 15 GHz and 22 GHz to be 
Ti5 0.2 and T22 ~ 0.05 respectively for component "A". 

No significant linear polarization is detected in the core re- 
gion at any of our frequencies with an approximate upper hmit 
of 1% at 22 GHz at the location of the strong, positive CP sig- 
nal analyzed here. The lack of significant hnear polarization 
is most likely explained by external depolarization in a thermal 
Faraday screen (Wardle 1971) which will have a negligible ef- 
fect on the circular polarization (Jones & O'DeU 1977; Homan 
etal. 2001). 

' Our procedure differs from a traditional CLEAN-type algorithm in that we improve the calibration with each cycle. We found a loop gain of unity to work well in 
these cases and tests with a loop gain of one-half showed no significant difference in the results. 



Relative R.A. (milli— arcsec) 

Fig. 1.— VLBA Stokes / images of 3C84 at 15 GHz. The base contour 
is 10 mJy/beam, and the contours increase in steps of \/2. The peaks inten- 
sity in the image is 2.79 Jy/beam. The FWHM dimensions of the naturally 
weighted beam (0.57 X 0.80 masxmas) is depicted by a cross figure in the 
lower left-hand corner of each panel. 

declining brightness profile, so Zero-V self-calibration will not 
gready distort the real Stokes V (Homan & Wardle 1999); this 
fact makes the reconstruction process simple and robust. 

Figure 2 shows the core region of 3C 84, the distorted CP- 
signal from Zem-V self-calibration, and the reconstruction of 
the true CP distribution. The reconstruction was done by as- 
suming the peak CP in the distorted map was real, subtracting 
this CP peak from the original data (prior to self-cahbration as- 
suming V = 0), calibrating the new subtracted data-set assuming 
y = 0, and examining the resulting residual map. After one it- 
eration, the residual V map still showed signs of real CP which 
had been distorted by assuming V = 0. (We expect this because 
the original subtraction could not have been large enough. The 
peak that was subtracted had already been suppressed by some 
factor in the initial calibration assuming V = 0.) So this pro- 
cess was repeated through a number of iterations until the final 
residual map was flat. We then restored all the subtracted com- 
ponents to the calibrated data. 

We ran a large number of tests on this CLEAN-type algo- 
rithm^ using simulated data, and we found the reconstruction to 
be robust, particularly for the bright positive components in fig- 
ure 2. The negative features are real; however, our tests on sim- 
ulated data showed that the reconstruction algorithm typically 
reduces their amplitude by ~ 10% at 15 GHz and ~ 30-40% at 
22 GHz and may introduce smaU shifts in the their peak posi- 
tion of ~ 0. 1 mas. These effects are both due to the weaker na- 
ture of negatively polarized components and their closer prox- 
imity to the peak of the total intensity structure where the most 
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Fig. 1 . — Core region of 3C 84 at 22 GHz (a) and 15 GHz (b). Left: Stokes / with contours beginning at 10 mjy/beam and increasing in steps of \/2. Center: 
Stokes V, calibrated assuming V = 0, creating some distortion to the true V distribution. Right: Stokes V, reconstruction of the true distribution using the algorithm 
described in §3.1. The V peaks in the reconstructed images are 12.7 mJy/beam at 22 GHz and 33.8 mJy/beam at 15 GHz. For aU Stokes-V images the contours begin 
at 2 mJy/beam and increase in steps of \/2. A single Stokes / contour is placed around the V images to show registration. The IJWHM dimensions of the naturally 
weighted beams ((0.42 x 0.59 masxmas at 22 GHz; 0.57 x 0.80 masxmas at 15 GHz) are depicted by a cross figure in the lower left-hand comer of each panel. 
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Fig. 3. — Circular polarization of core region of 3C 84 with matched beams 
at 15 and 22 GHz. The position,s of three circularly polarized components are 
marked in the 22 GHz image. These same positions are indicated in the 15 GHz 
image to show registration. A single Stokes / contour is placed around the V 
images to show registration. The common FWHM restoring beam (0.47 x 0.63 
masxmas) is depicted by a cross figure in the lower left-hand comer of the im- 
age. 

4. DISCUSSION AND CONCLUSIONS 

We first consider the spot of positive circular polarization 
which is extremely strong at -1-3% of the local Stokes / at 15 
GHz. If the circular polarization is intrinsic to the synchrotron 
radiation itself, its fractional level varies as oc v~^^-^ for opti- 
cally thin emission. This is flatter than our observed spectrum 
but differs from it by only l.Scr. The full expression is (Jones 
& O'Dell 1977; Wardle & Homan 2003) 



0.5 "uu 



(1) 



where is a constant « 2 for a < -1 (Jones & O'Dell 1977) 
and ub^^ = 2.85™ MHz is the gyro-frequency for Bj™ in 
Gauss. BuuCosO is the component of uniform, unidirectional 
magnetic field along the line of sight responsible for generating 
the circular polarization. B™* is the RMS averaged component 
of field in the plane of the sky and includes contributions not 
only from the transverse part of B„j, but also from disordered 
field or ordered loops of field which will contribute to Stokes / 
but not V. 

If the circular polarization is entirely due to the intrin- 
sic mechanism, we can estimate the minimum magnetic field 
strength required by the very high levels of circular polarization 
detected at 15 GHz. Walker, Romney, & Benson (1994) find the 
jet of 3C 84 to be only mildly relativistic and oriented at a large 
angle to the line of sight, 6 = 30° -55°. We take the Doppler 
factor to be unity and 9 = 30° in these calculations. The largest 
possible circular polarization will be produced by the intrin- 
sic mechanism when the magnetic field is completely uniform 
and uni-directional. Under these conditions, B™^ = sin 9, 
and we find that a uniform field strength of B„„ ~ 1 Gauss is 
necessary to produce the observed 3.2% circular polarization 
at = 15 GHz. Any sign reversals or disorder in the magnetic 
field will increase the required magnetic field strength. 

Given the above estimate, it is possible that a very strong 
(~ 1 G), uniform magnetic field along the jet axis is responsi- 
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ble for generating the observed positive circular polarization via 
the intrinsic mechanism alone. However, the negative compo- 
nents of circular polarization cannot be easily explained by the 
intrinsic mechanism. They are located upstream where the jet 
is optically thick. Pacholczyk & Swihart (IQ?!)*" showed that 
intrinsic circular polarization will change sign when the optical 
depth is significantly greater than unity, but their results also 
show that the spectral index of the fractional polarization, nic, 
can only be negative after the sign change. This result conflicts 
with our observations which show a positive spectral index for 
nic for both negatively polarized components. 

An alternative to a pure intrinsic circular polarization model 
is to have a significant (and perhaps dominant) contribution 
from the Faraday conversion effect. High levels of circular po- 
larization can be easily generated via the Faraday conversion 
process which converts linear polarization to circular (Jones 
1988). Faraday conversion operates on Stokes U relative to 
the local magnetic field orientation (while synchrotron radia- 
tion from the same field generates only Stokes Q) and therefore 
requires either some change in field orientation along the Une of 
sight or some internal Faraday rotation to drive the process. If 
the conversion process is driven by Faraday rotation, some uni- 
form, unidirectional magnetic field, Z?„„, is required to generate 
internal Faraday rotation, but here we can have <C 5™. 

In the optically thin regime, the conversion process has a 
steep spectrum, oc z^~^, or even steeper if it is driven by 
Faraday rotation (Wardle & Homan 2003). For the positively 
polarized, optically thin component, our observed spectrum is 
flatter than this, jy-c^io ^^ ^jj^j yyould seem to reject Faraday 
conversion as the dominant process. However, at ris ~ 0.2, we 
are not completely optically thin, and the predicted spectrum 
flattens as we approach r = 1 (Jones & O'Dell 1977). This is 
due to optical depth and saturation in the conversion process 
as it cycles into a regime where it generates the opposite sign 
of circular polarization, depolarizing the radiation and flatten- 
ing the spectrum. This saturation most easily occurs when the 
conversion process is driven by internal Faraday rotation, as the 
Faraday rotation will also saturate and begin to depolarize the 
linear polarization, reducing the amount available for conver- 
sion and also flattening the spectrum. 

Detailed modeling of this saturation effect is beyond the 
scope of this paper; however, it is possible to explain not only 
the spectrum of the positively polarized component but also the 
change in sign and the spectrum of the negatively polarized 
components via this mechanism. Of course, the sign change 



of polarization from negative to positive may not be due to sat- 
uration but rather due to changes in magnetic field order. If the 
Faraday conversion is driven by Faraday rotation, it is possible 
to flip the net field polarity and generate the opposite sign of 
circular polarization. If the Faraday conversion is not driven by 
Faraday rotation but rather results from a changing field order 
along the line of sight, as might be the case for a helical field, 
a change in that field order, such as the pitch angle of the helix, 
could also change the sign of the resulting circular polarization. 

Another interesting aspect of the circular polarization dis- 
tribution is its apparent break-up into three distinct circularly 
polarized components while the total intensity (Stokes-/) dis- 
tribution appears smooth and continuous throughout the region. 
The sign change of circular polarization between components A 
and B will naturally make those regions appear separated, since 
Stokes-y must pass through zero between them. The discrete 
nature of component C is more difficult to understand. We note 
that it is closer to the base of the jet and offset from its center, 
so its discrete appearance may be due to a fortuitous location 
where the optical depth is not too large and other parameters 
are suitable for producing significant CP. The magnetic field 
order in the jet and the presence (or absence) of shocks in the 
magnetic field may also create enhanced regions of circular po- 
larization; however, without linear polarization information it 
is impossible to evaluate those possibilities. 

Finally, we note that 3C 84, at a redshift of z = 0.017, shows 
the strongest fractional circular polarization yet detected by 
the VLBA. Much more powerful quasars rarely show local 
CP stronger than 0.3% (Homan & Wardle 1999; Homan et al. 
2001); however, these sources are much farther away and at 
such distances 3C 84 would have only 0.1 -0.2% local circu- 
lar polarization due to blending with the core. In this light, it is 
interesting that the intra-day variable source PKS 1519-273 ex- 
hibits CP at the 2-4% level in a scintillating component 15-35 
/zas in size (Macquart et al. 2000). This corresponds to a region 
of 0.05-0.30 pc in size (for z > 0.2, Stickel, Fried, & Kuehr 
1993), comparable to the linear resolution we achieve in these 
observations, ~ 0.15 pc. The coincidence of high levels of CP 
with small spatial scales may be due to to the increased relative 
importance of uni-directional magnetic field (£«„) in the jet at 
smaU radii (Wardle & Homan 2003). 
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* The plots given by Jones & O'Dell (1977) consider the combined effects of intrinsic circular polarization and Faraday conversion and here we are interested in the 
intrinsic mechanism alone. We verified the Pacholczyk & Swihart (1971) result by direct numerical integration of the Jones and O'Dell radiative transfer formulas. 



